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Preface 


The  two-dimensional  numerical  modeling  of  the  flow  conditions  in  the 
proposed  Kentucky  Lock  was  performed  for  the  U.S.  Army  Engineer  District, 
Nashville  (LRN).  This  study  was  authorized  by  the  U.S.  Army  Engineer 
Division,  Great  Lakes  and  Ohio  River,  on  13  January  1997.  Mr.  Donald  Getty, 
LRN,  directed  this  study. 

This  work  was  conducted  in  the  Coastal  and  Hydraulics  Laboratory  (CHL) 
of  the  U.S.  Army  Engineer  Waterways  Experiment  Station  (WES)  during  the 
period  March  1997  to  October  1997  under  the  direction  of  Dr.  J.  R.  Houston, 
Director,  CHL;  Mr.  Charles  C.  Calhoun,  Assistant  Director,  CHL;  and  Dr.  P.  G. 
Combs,  Chief,  Rivers  and  Structures  Division,  CHL. 

Extension  of  the  HTVEL2D  (two-dimensional  hydrodynamic  numerical) 
model  was  completed  by  Dr.  R.  L.  Stockstill,  Spillways  and  Channels  Branch, 
Rivers  and  Structures  Division,  CHL.  Simulation  runs  and  analyses  of  results 
were  conducted  by  Dr.  J.  E.  Hite,  Jr.,  Leader  of  the  Locks  and  Conduits  Group, 
Rivers  and  Structures  Division,  and  Dr.  Stockstill  under  the  supervision  of 
Mr.  B.  P.  Fletcher,  Chief,  Spillways  and  Channels  Branch.  Mr.  S.  Cornell, 
Spillways  and  Channels  Branch,  also  assisted  with  the  simulations  and 
postprocessing  of  the  model  results.  Dr.  R.  C.  Berger,  Estuaries  and 
Hydroscience  Division,  CHL,  provided  technical  assistance  and  peer  review. 

The  report  was  written  by  Drs.  Hite  and  Stockstill. 

At  the  time  of  publication  of  this  report.  Director  of  WES  was  Dr.  Robert  W. 
Whalin.  Commander  was  COL  Robin  R.  Cababa,  CE. 


The  contents  of  this  report  are  not  to  be  used  for  advertising ;  publication  or 
promotional  purposes.  Citation  of  trade  names  does  not  constitute  an 
official  endorsement  or  approval  of  the  use  of  such  commercial  products. 


Conversion  Factors,  Non-SI  to 
SI  Units  of  Measurement 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI  units 
as  follows: 


|  Multiply 

By 

To  Obtain  | 

|  cubic  feet 

0.02831685 

cubic  meters  1 

I  cubic  yards 

0.7645549 

cubic  meters  I 

feet 

0.3048 

meters  8 

miles  (U.S.  statute) 

1.609347 

kilometers  fl 

tons  (force) 

8.896443 

kiionewtons  I 
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Introduction 


Background 

The  existing  Kentucky  Lock  is  located  on  the  Tennessee  River  approximately 
20  miles1  southeast  of  Paducah,  KY  (Figure  1).  The  project  consists  of  a  gated 
spillway  to  regulate  riverflows,  a  powerhouse  for  hydroelectric  power 
generation,  and  a  600-ft-long  navigation  lock  for  moving  industrial  tow  traffic 
and  recreational  boats  through  the  project  The  existing  lock  is  operating  at 
capacity,  and  an  additional  1,200-ft-long  by  110-ft-wide  lock  is  projected  to  be 
necessary  to  satisfy  future  capacity  requirements.  The  new  lock  will  be  located 
landward  of  the  existing  lock  with  the  upstream  pintles  (cross-stream  axis  of  die 
miter  gates)  located  just  over  100  ft  downstream  from  the  upstream  pintles  of  the 
existing  lock. 

The  new  lock  features  a  through-the-sill  intake  that  carries  flow  to  a  multiport 
filling  and  emptying  system.  One  proposed  lock  discharge  plan  uses  an 
interlaced  lateral  system  located  downstream  of  the  lower  miter  gate  pintle  as 
shown  in  Figure  2.  Another  alternative  being  investigated  for  the  discharge 
system  is  a  landside  channel  that  discharges  downstream  of  the  lower  approach 
guide  wall  as  shown  in  Figure  3.  A  final  decision  on  the  discharge  system  will 
be  madft  after  all  alternatives  are  evaluated. 


Purpose  and  Scope 

The  initial  investigation  was  performed  to  evaluate  the  flow  conditions  in  die 
lower  lock  approach  for  different  excavation  plans  in  the  approach  channel  with 
the  interlaced  lateral  discharge  system.  These  flow  conditions  must  be  known  to 
determine  the  effect  these  flows  have  on  tows  in  the  lower  approach  area. 
Adverse  flow  conditions  (large  streamwise  and  cross-stream  water-surface 
gradients)  in  the  lower  approach  may  prohibit  tows  from  mooring  in  this  area 
during  lock  discharges.  The  flow  conditions  selected  for  evaluation  were  a 


1  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI  units  is  found  on  page  vi. 
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headwater  el1  of  357,  a  tail  water  el  of  304.2,  and  emptying  valve  opening  times  of 
1.5  min  and  1 1.7  min.  This  flow  is  considered  the  50  percent  duration  condition, 
and  the  valve  speeds  represent  fast  and  slow  valve  operations,  respectively. 

Two  excavation  plans  were  evaluated  with  the  interlaced  lateral  discharge 
design:  the  minimum  excavation  plan  and  the  moderate  excavation  plan.  The 
invert  elevation  of  the  lower  approach  channel  with  the  minimum  excavation  plan 
was  289,  and  with  the  moderate  excavation  plan  was  284. 

A  subsequent  investigation  of  a  landside  discharge  channel  was  performed 
after  evaluation  of  the  interlaced  lateral  discharge  system.  The  U.S.  Army 
Engineer  District,  Nashville,  requested  that  the  U.S.  Army  Engineer  Waterways 
Experiment  Station  (WES)  evaluate  the  flow  conditions  in  the  lower  lock 
approach  with  a  headwater  el  of  359  and  tailwater  el  of  302  and  3 14  with  an 
emptying  valve  operation  of  13  min.  Also,  flow  conditions  were  evaluated  with 
a  headwater  el  of  359,  a  tailwater  el  of  302,  and  an  emptying  valve  operation  of 
6  min. 

The  landside  and  interlaced  lateral  discharge  alternatives  were  evaluated  by 
comparing  water-surface  differentials  at  selected  locations.  Hie  hawser  forces  a 
tow  and  barge  arrangement  will  experience  are  directly  related  to  the  water- 
surface  slope  on  which  die  vessel  rests.  The  particular  locations  for  the 
interlaced  lateral  discharge  alternative  were  near  the  bow  and  stem  of  a  3  x  3 
barge  arrangement  and  a  3  x  5  barge  arrangement  both  having  their  upstream  end 
moored  at  sta  17+00.  Flow  conditions  resulting  from  the  landside  discharge 
alternative  were  also  evaluated  at  sta  15+00  and  27+00,  which  would  be 
representative  of  a  3  x  6  barge  arrangement  with  the  upstream  end  moored  at 
sta  15+00. 


Approach 

The  two-dimensional  (2D),  depth-averaged  flow  model,  HTVEL2D,  was  used 
to  model  the  unsteady  velocities  and  water-surface  elevations  in  die  lower  lock 
approach  resulting  from  lock  emptying  operations.  The  HTVEL2D  model  was 
chosen  for  this  study  because  it  is  designed  to  provide  numerically  stable 
solutions  for  advection-dominated  flow  containing  large  gradients  in  the  flow 
variables.  Large  gradients  in  depths  and  velocities  are  present  in  the  vicinity  of 
lock  outlets  during  emptying  operations.  The  flow  conditions  in  these  areas  can 
vary  from  no  flow  to  peak  discharges  of  about  22,000  cfs  in  less  than  2  min  The 
HTVEL2D  code  was  modified  to  allow  specification  of  time-dependent  inflow 
boundary  conditions. 

A  plan  view  of  the  geometry  and  topography  in  the  lower  lock  approach  with 
the  interlaced  lateral  discharge  alternative  and  the  minimum  excavation  plan  is 


1  All  elevations  (el)  and  stages  cited  herein  are  in  feet  referred  to  the  National  Geodetic  Vertical 
Datum  (NGVD). 
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shown  in  Figure  4.  The  mesh  used  for  the  hydrodynamic  computations  with  the 
initial  simulations  is  shown  in  Figure  5.  The  computational  mesh  used  for  the 
simulations  with  the  landside  discharge  channel  is  shown  in  Figure  6. 
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2  Description  of  Model 


Governing  Equations 

Fluid  motion  is  modeled  using  the  2D  unsteady  shallow-water  equations.  The 
shallow-water  (or  long-wave)  equations  are  a  result  of  die  vertical  integration  of 
the  equations  of  mass  and  momentum  conservation  for  incompressible  flow 
under  the  hydrostatic  pressure  assumption.  This  assumption  implies  that  vertical 
accelerations  are  negligible  compared  with  the  horizontal  accelerations  and  the 
acceleration  due  to  gravity.  The  vertical  accelerations  are  small  when  the 
characteristic  wavelength  is  long  relative  to  the  depth,  which  is  why  these 
equations  are  referred  to  as  long-wave  or  shallow-water  equations. 

The  dependent  variables  of  the  fluid  motion  are  defined  by  die  flow  depth  k, 
the  x-direction  component  of  unit  discharge  p,  and  the  y-direction  component  of 
unit  discharge  q.  These  variables  are  functions  of  the  independent  variables  x 
and  y,  the  two  space  directions,  and  time  t.  If  the  free-surface  stresses  are 
neglected,  the  shallow-water  equations  are  given  as  (Abbott  1979): 
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G  = 


(4) 


and 


H  = 


,  ^0  2  P^Ip*  + 

^  +  Tl^gi^ - 3_ 

ox 


Coh 


7/3 


,  2  q>Ip*  +  ?2 

«*-z-  +  *  g'  '  7n 

5y  c^7* 


(5) 


where 


g  =  acceleration  due  to  gravity 
p  =  fluid  density 
Zg  =  channel  bed  elevation 
n  =  Manning's  roughness  coefficient 


C0  =  dimensional  constant  (C0  =  1  for  SI  units  and  2.208  for 
non-SI  units) 

Oyp  Oyy,  =  Reynolds  stresses  due  to  turbulence  where  the  first  subscript 
indicates  the  direction  and  the  second  indicates  the  face  on 
which  the  stress  acts 

The  Reynolds  stresses  are  determined  using  the  Boussinesq  approach  relating 
stress  to  the  gradient  in  the  mean  currents: 


(6) 


yy 


(7) 


and 
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(8) 


yx 


=  PVr 


du  dv 
dy  dx 


where 

vf  =  kinematic  eddy  viscosity,  which  varies  spatially 

u  =  p/h,  the  depth-averaged  x-direction  component  of  velocity 

v  =  q/h,  die  depth-averaged  y-direction  component  of  velocity 

Values  of  the  eddy  viscosity  are  determined  empirically  as  a  function  of  the  local 
flow  variables  as  (Rodi  1980;  Chapman  and  Kuo  1985): 

v,  =  «v/8i£  h5'6  \V\  (9) 

co 


where 

C  =  coefficient  that  varies  between  0.1  and  1.0 
\V\  =  velocity  vector  magnitude  =  (h2  +  v2)1® 


Extensions  to  Include  Unsteady  Inflow 

HTVEL2D  version  1.07  (Stockstill  and  Berger  1994)  was  designed  to  serve  as 
a  tool  to  provide  steady-state  solutions  of  flow  fields.  Modeling  flow  conditions 
in  the  lower  lock  approach  during  lock  emptying  required  the  ability  to  simulate 
the  time-dependent  flow  rates  that  represent  the  lock  emptying  hydrograph  in  the 
computational  model.  The  time-dependent  discharge  from  the  interlaced  laterals 
was  the  inflow  used  to  drive  the  hydrodynamic  model.  Two  methods  were  used 
to  simulate  the  lock  discharge.  In  the  first  method,  the  total  discharge  was 
applied  as  flux  through  an  inflow  boundary  located  just  downstream  of  the 
discharge  laterals.  This  meant  that  values  of  p  and  q  at  this  inflow  boundary 
were  specified  as  functions  of  time.  Since  significant  vertical  accelerations  are 
generated  in  the  real  system  in  the  vicinity  of  the  discharge  laterals,  there  was 
concern  that  the  horizontal  momentum  resulting  from  these  specified  boundary 
conditions  was  unrealistic.  Therefore,  a  second  method  was  developed  to 
simulate  the  lock  discharge. 

In  this  method  the  mass  flow  rate  resulting  from  the  lateral  discharge  outlet 
area  was  specified  as  a  point  source  of  mass.  The  conservation  of  mass  equation 
(first  row  vector  in  Equation  1)  was  extended  to  include  a  change  in  mass  term 
on  the  right-hand  side. 
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(10) 


dh  ^  dp  <  dq  _  . 
5/  Sr  dy 


where  i  is  a  time-dependent  scalar  value  having  units  of  length  per  time  Values 
of  z  are  determined  such  that 

/  m  dAi  =  Q{i)  (11 

A 


where 

z'ft)  =  time  dependent  source  term 

At  =  plan  area  of  the  interlaced  lateral  structure 

Q( t)  =  lock  discharge  from  the  emptying  hydrograph 

This  method  of  representing  the  discharge  hydrograph  was  considered  an 
improvement  over  the  first  method  because  it  did  not  impose  horizontal 
momentum  at  the  inflow  boundary.  Numerical  model  results  revealed 
insignificant  differences  in  the  flow  variables  two  computational  elements  away 
from  the  discharge  laterals.  Peak  water-surface  elevations  computed  with  the 
point  source  method  slightly  lagged  in  time  behind  the  peaks  computed  when 
horizontal  momentum  was  applied  at  die  inflow  boundary.  The  two  methods  of 
mass  input  provided  essentially  the  same  results  at  locations  away  from  the 
discharge  laterals.  The  point  source  method  better  represented  die  real  system 
since  no  horizontal  momentum  was  assumed  and  currents  were  generated  by  the 
water-surface  gradients  across  the  area  representing  the  discharge  laterals  and  the 
remaining  flow  field. 


Discretization 

Because  there  are  no  general  solutions  of  the  continuous  functions  in  the 
governing  equations,  they  must  be  solved  by  numerical  methods.  Discrete  values 
of  die  unknown  variables  are  solved  using  a  Petrov-Galerkin  finite  element 
representation  of  the  equations.  Details  of  the  finite  element  formulation  are 
presented  in  Stockstill  and  Berger  (1994). 

A  finite  difference  expression  is  used  for  the  temporal  derivatives.  The 
general  expression  for  die  temporal  derivative  of  the  unknowns,  U},  is: 

=  —(uM  -  u.k)  +  —  g)  (uk  -  Uk~l)  (12) 

dt  At '  J  J  1  At  '  3  1  I 
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where 


a  =  temporal  differencing  weight 
At  =  time-step  size 
j  =  node  location 
k  -  time-step 

An  a  value  of  1  results  in  a  first-order  backward  differencing,  and  an  a  of  1.5 
results  in  a  second-order  backward  differencing  approximation  of  the  temporal 
derivative. 
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Model  Applications 


Application  of  the  modified  HTVEL2D  model  to  the  proposed  lower  approach 
required  the  construction  of  a  numerical  model  computational  mesh  to  represent 
the  design  geometry  and  the  specification  of  boundary  conditions  and  model 
parameters. 


Geometry  and  Computational  Mesh 

Interlaced  lateral  discharge  alternative  computational  mesh 

Initially,  the  lower  approach  length  to  the  new  lock  was  modeled  from  the 
upstream  end  of  the  interlaced  laterals  (approximately  sta  14+50)  to  130  ft 
downstream  from  the  end  of  the  proposed  guide  wall  (sta  31+00).  The  entire 
width  of  the  lower  approach  was  modeled  down  to  sta  25+50;  and  from  sta 
25+50  to  31+00,  a  700-ft  width  of  the  lower  approach  was  modeled.  Figure  4 
shows  details  of  the  geometry  of  the  lower  lock  approach.  This  initial  mesh  had 
564  nodes  and  609  elements.  Grid  resolution  was  more  refined  at  alignment  and 
grade  changes. 

Examination  of  model  results  with  the  initial  mesh  led  to  the  conclusion  that 
the  downstream  boundary  condition  had  significant  influence  on  the  velocity 
distribution  on  the  left  bank  at  the  downstream  end  of  the  approach  channel. 
Therefore,  the  model  limits  were  extended  significantly  to  prevent  the  outflow 
boundary  conditions  from  influencing  tine  flow  patterns  at  the  end  of  the 
approach  channel.  The  larger  mesh,  shown  in  Figure  5,  included  a  4,400-ft 
length  of  the  Tennessee  River.  The  entire  width  of  the  river  was  represented 
from  the  tailrace  of  the  powerhouse  and  spillway  to  4,700  ft  downstream  from 
the  outlet  This  mesh  had  1012  nodes  and  1529  elements.  Resolution  was 
increased  where  large  gradients  of  the  flow  variables  were  expected. 


Landside  discharge  channel  alternative  computational  mesh 

This  alternative  required  a  large  computational  domain  so  that  the  flow 
patterns  were  properly  modeled  near  the  vicinity  of  the  discharge  channel  outlet 
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This  computational  mesh  had  1 140  nodes  and  1851  elements  and  was  finely 
resolved  at  the  discharge  channel  outlet  to  the  river  as  shown  in  Figure  6. 


Boundary  Conditions 

Interfaced  lateral  alternative  boundary  conditions 

This  study  included  simulation  of  two  discharge  hydrographs  furnished  by  the 
Nashville  District  that  represented  emptying  valve  openings  of  1.5  min  and 
1 1.7  min.  The  hydrographs  for  these  valve  operations  are  shown  in  Figure  7. 
Tables  1  and  2  provide  the  discharge  and  time  data  for  these  hydrographs. 


Landside  discharge  channel  alternative  boundary  conditions 

Inflow  boundary  conditions  for  this  study  were  time-histories  of  the  flow  rate 
at  sta  26+20  of  the  landside  discharge  channel.  The  two  hydrographs  shown  in 
Figure  8  represented  an  emptying  valve  time  of  1-5  min  with  a  headwater  el  of 
359  and  tailwater  el  of  302  and  314.  Tables  3  and  4  provide  the  discharge  and 
time,  data  for  these  hydrographs.  The  third  condition  evaluated  with  the  landside 
discharge  channel  alternative  was  a  headwater  el  of  359,  a  tailwater  el  of  302, 
and  a  6-min  emptying  valve.  The  hydrograph  for  these  conditions  is  shown  in 
Figure  9,  and  the  discharge  and  time  data  are  provided  in  Table  5.  As  with  the 
interlaced  lateral  design,  the  hydrographs  used  as  inflow  boundary  conditions  for 
the  landside  discharge  channel  were  furnished  by  the  Nashville  DistricL 


Modei  Parameters 


Model  and  flow  parameters  used  in  the  simulations  are  provided  in  the 
following  tabulation. 


|  Model  Condition 

Value  E 

a 

1-25  8 

3 

025  1 

2.208  1 

c 

as  | 

n 

0.02  | 

B  ......  .... 

32.2  ft/sec2 

Here,  {5,  which  can  range  in  value  from  0.0  to  0.5,  is  a  weighting  term  used  in 
the  test  function  of  the  finite  element  formulation  of  the  governing  equations. 
Details  of  the  model’s  finite  element  formulation  are  provided  in  Stockstill  and 
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Berger  (1994).  P  can  be  thought  of  as  an  upwinding  parameter  that  provides  the 
numerical  stability  needed  for  modeling  advection-dominated  flows  that  may 
have  shocks  such  as  hydraulic  jumps.  Actually,  HTVEL2D  allows  the  input  of 
two  values  of  P  (Berger,  Stockstill,  and  Ott  1995).  Small  values  (i.e.,  0.1  to  0.2) 
are  more  precise  and  have  been  successfully  applied  to  regions  of  the  flow  field 
having  smooth  solutions.  Therefore,  the  larger  value  (i.e.,  0-5)  is  automatically 
applied  only  to  the  domain’s  roughest  regions,  which  are  generally  shocks  such 
as  hydraulic  jumps  (Berger  and  Stockstill  1995;  Stockstill,  Berger,  and  Nece 
1997).  The  computed  flow  field  was  not  expected  to  contain  any  hydraulic 
jumps;  therefore  only  one  value  of  P  was  used,  which  was  chosen  as  an  average 
of  values  used  for  smooth  and  rough  solutions. 

Berger  (1993),  while  investigating  the  modeling  of  dam  breaks,  found  that  a 
value  of  1.25  for  the  temporal  derivative  weighting  coefficient,  a,  provides 
accurate  timing  of  problems  involving  rapidly  varying  flow.  Choice  of  the 
turbulent  eddy  viscosity  coefficient,  C,  was  simply  an  average  of  the  range  over 
which  this  coefficient  has  been  known  to  vary  (0.1  to  1.0) .  Solutions  are 
relatively  insensitive  to  the  selection  of  this  parameter  since  the  Boussinesq 
stresses  describe  momentum  dispersion,  which  has  insignificant  effects  on  the 
computed  flow  depths.  Likewise  the  Manning  coefficient  n  of  0.02  probably 
describes  a  channel  that  is  slightly  smoother  than  fire  rock-cut  channel  proposed 
for  die  Kentucky  Lock,  but  water-surface  differentials  within  the  channel 
generated  by  the  rapidly  varying  flow  issuing  from  the  lock  outlets  are 
dominated  by  wave  characteristics  and  not  bed  drag. 
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4  Results 


The  flow  conditions  in  the  lower  approach  were  evaluated  by  comparing  the 
water-surface  elevations  at  two  selected  locations.  For  example,  wife  fee 
interlaced  lateral  discharge  alternative,  fee  upstream  location  was  147.5  ft 
downstream  from  the  downstream  end  of  fee  interlaced  lateral  (sta  17+00)  and 
55  ft  out  from  fee  guide  wall  (Figure  4).  One  of  fee  downstream  locations  was 
600  ft  from  fee  upstream  location  (sta  23+00)  and  also  55  ft  out  from  fee  guide 
wall  In  this  example,  these  water  surfaces  would  represent  fee  water-surface 
slope  in  fee  lower  approach  to  which  a  3  x  3  barge  arrangement  moored  142_5  ft 
downstream  from  fee  interlaced  laterals  would  be  exposed.  A  tow  and  barges 
moored  in  this  area  and  exposed  to  extreme  water-surface  slopes  and  high- 
velocity  flows  will  experience  large  hawser  forces. 


Interlaced  Lateral  Discharge  Alternative 

Simulations  wife  both  excavation  plans  were  conducted  using  a  time-step  of 
4  sec.  The  total  time  modeled  for  each  of  the  valve  schedules  was  feat  given  by 
fee  lock  discharge  hydrographs  feat  served  as  inflow  boundary  conditions. 
Longer  simulation  times  were  not  run  because  inflow  boundary  conditions  were 
unknown  beyond  fee  times  given  by  fee  hydrographs. 


Minimum  excavation  plan 

1.5-min  valve,  sta  17+00  and  sta  23+00.  The  initial  simulation  was 
performed  wife  fee  minimum  excavation  plan,  headwater  el  of  357,  tailwater  el  of 
304.2,  and  an  emptying  valve  opening  time  of  1.5  min.  These  headwater  and 
tailwater  conditions  were  performed  for  all  simulations  wife  fee  interlaced  lateral 
discharge  alternative.  In  future  descriptions  of  flows  wife  the  interlaced  laterals 
these  conditions  are  assumed-  A  time-history  of  fee  water-surface  elevation  at 
sta  17+00  is  shown  in  Figure  10.  The  maximum  water  surface  was  305.7  (which 
was  13  ft  higher  than  the  tailwater)  and  occurred  1 12  sec  after  fee  valve  began 
opening.  The  minimum  water-surface  elevation  at  sta  17+00, 55  ft  out  from  the 
guide  wall,  was  302.7  and  occurred  8.8  min  after  fee  emptying  valve  began 
opening. 
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The  water-surface  elevations  increased  in  the  lower  approach  as  the  positive 
wave  produced  by  the  lock  discharge  traveled  downstream.  For  this  discussion, 
“positive”  was  used  to  designate  a  wave  surface  higher  than  the  initial  steady- 
state  surface.  If  the  wave  surface  was  lower  than  the  steady-state  surface,  the 
wave  was  described  as  negative.  A  wave  reflected  at  the  upstream  end  (lock 
miter  gates)  by  doubling  in  magnitude  when  it  intercepted  the  wall.  When  a  wave 
reached  the  abrupt  expansion  in  the  lower  approach,  it  was  “pinned”  in  height 
and  therefore  reflected  as  a  negative  wave.  This  produced  a  standing  wave  (the 
entire  approach  length  rocked  up  and  down  in  phase).  The  seiche  frequency  was 
determined  by  the  geometry  of  the  approach  and  the  depth.  Several  projects  (e.g.. 
Shows  and  Franco  1981  and  Oswalt,  Abies,  and  Murphy  1972)  experience 
oscillations  within  approach  canals  when  the  lock  inrakes  or  outlets  are  located  in 
or  near  the  channel 

A  time-history  of  the  water-surface  elevation  at  the  downstream  location 
(sta  23+00, 55  ft  out  from  the  guide  wall)  is  also  shown  in  Figure  10.  The 
maximum  water  surface  is  slightly  lower  (305.3)  than  that  observed  at  sta  17+00 
and  the  minimum  water-surface  elevation  is  slightly  higher  (3033).  A  time- 
history  of  the  water-surface  differential  between  sta  17+00  and  sta  23+00, 55  ft 
from  the  guide  waft,  is  shown  in  Figure  1 1.  The  downstream  differential  (which 
is  considered  a  positive  water-surface  slope)  results  from  water  surfaces  at 
sta  17+00  being  higher  than  those  at  sta  23+00  at  the  same  time,  and  the 
upstream  differential  (a  negative  water-surface  slope)  results  from  water  surfaces 
at  sta  17+00  being  lower  than  those  at  sta  23+00  at  the  same  time.  A  positive 
water-surface  slope  would  move  an  unmoored  tow  and  barges  downstream,  and  a 
negative  water-surface  slope  would  result  in  an  upstream  drift  of  a  tow  and 
barges.  The  maximum  positive  slope  was  0.63  ft  in  600  ft,  and  the  maximum 
negative  slope  was  0.74  ft  in  600  ft  Neglecting  the  forces  due  to  drag  and 
inertia,  assuming  that  the  barge  train  acts  as  a  single  rigid  vessel  and  neglecting 
the  effect  of  the  vessel  blockage  area  of  the  approach  channel  the  force  required 
to  hold  a  vessel  in  place  is  a  function  of  only  the  water-surface  slope.  Using 
standard  barge  dimensions  of  195  ft  by  35  ft,  the  longitudinal  hawser  force 
computed  for  a  3  x  3  barge  arrangement  with  a  9-ft  draft  and  a  water-surface 
slope  of  -0.74  ft  in  600  ft  is  213  tons.  Due  to  the  negative  slope,  a  tow  and 
barge  arrangement  would  tend  to  move  toward  the  downstream  miter  gates. 
Realizing  that  these  assumptions  are  not  entirely  correct,  one  can  still  obtain  a 
relative  sense  of  forces  in  the  lower  lock  approach  due  to  water-surface  slopes. 

A  plan  view  of  velocity  vectors  and  water-surface  contours  is  shown  in 
Figure  12.  These  are  the  flow  conditions  that  occur  92, 180,  and  360  sec  after  the 
emptying  valve  began  opening.  The  highest  velocity  computed  for  sta  17+00, 

55  ft  from  the  guide  wall  occurred  at  204  sec  and  was  5  ft/sec. 

13-nun  valve,  sta  17+00  and  27+00.  Additional  time-histories  were  plotted 
to  evaluate  the  water-surface  slopes  between  sta  17+00  and  27+00.  These 
conditions  are  applicable  to  a  3  x  5  barge  arrangement  in  the  lower  approach.  The 
maximum  water-surface  differential  was  found  to  be  -1.19  ft  (Figure  1 1),  which 
would  result  in  a  longitudinal  hawser  force  of  34.2  tons  directed  toward  the  miter 
gates,  for  the  assumptions  previously  mentioned. 
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11.7-min  valve,  sta  17+00  and  23+00-  The  next  simulation  was  performed 
with  minimum  excavation  and  the  slower,  1 1 .7-min  valve  opening.  Time- 
histories  of  water-surface  elevation  with  the  11.7-min  valve  at  sta  17+00  and 
23+00, 55  ft  from  the  guide  wall,  are  shown  in  Figure  10.  The  maximum  water- 
surface  elevation  was  304.6  and  the  minimum  water-surface  elevation  was  303.8 
at  sta  17+00  as  shown  in  Figure  10.  The  water-surface  elevations  at  sta  23+00 
are  very  similar  to  those  at  sta  17+00  as  Figure  10  illustrates.  These  results 
indicate  that  with  the  slower  valve  opening,  the  water  surface  in  the  lower 
approach  varies  only  slightly  from  the  tailwater  elevation.  The  maximum  water- 
surface  differential  with  the  long  valve  time  (1 1.7  min)  was  0.10  ft  as  shown  in 
Figure  11.  If  tins  differential  is  converted  to  a  longitudinal  hawser  force  as 
described  previously,  the  computed  hawser  force  for  a  3  *  3  barge  arrangement 
with  the  upstream  end  moored  at  sta  17+00  is  2.9  tons.  Flow  conditions  in  the 
lower  lock  approach  were  improved  with  the  slower  valve  speed;  however,  lock 
emptying  times  will  increase  with  the  slower  valve. 

11.7-min  valve,  sta  17+00  and  27+00.  Additional  time-histories  were 
obtained  with  the  1 1 .7-min  valve  speed  to  evaluate  the  water-surface  slopes 
between  sta  17+00  and  27+00.  The  maximum  water-surface  differential  was 
found  to  be  -0.16  ft  (Figure  1 1),  which  would  result  in  a  longitudinal  hawser 
force  of  4.6  tons  with  the  assumptions  stated  previously. 

Comparison  of  results  with  minimum  excavation.  A  comparison  of  the 
water-surface  differentials  at  sta  17+00  and  23+00, 55  ft  from  the  guide  wall, 
between  the  1.5-  and  11.7-min  valve  times  is  shown  in  Figure  1 1.  The  results 
clearly  indicate  that  a  slower  valve  time  causes  much  less  water-surface  slope  in 
the  lower  approach  channel  and  would  be  more  favorable  for  tows  and  barges 
trying  to  moor  in  this  area.  This  is  also  apparent  when  comparing  the  results 
between  sta  17+00  and  27+00  shown  in  Figure  1 1. 


Moderate  excavation  plan 

1-5-min  valve,  sta  17+00  and  23+00.  The  next  simulation  was  performed 
with  the  moderate  excavation  plan  and  a  1-5-min  opening  schedule  for  the 
emptying  valves.  The  moderate  excavation  plan,  shown  in  Figure  13,  consisted 
of  excavating  approximately  57,000  yd3  in  an  area  1,500  ft  long  by  200  ft  wide 
beginning  just  downstream  from  the  interlaced  laterals.  The  bed  elevation  in  this 
area  was  generally  5  ft  lower  than  that  of  the  minimum  excavation  plan. 

Time-histories  of  the  water-surface  elevations  at  sta  17+00  and  23+00  with 
the  1.5-min  valve  and  the  moderate  excavation  plan  are  shown  in  Figure  14. 

The  maximum  water-surface  at  sta  17+00  occurred  104  sec  after  the  valve  began 
opening  and  reached  el  305.6,  which  is  1.4  ft  higher  titan  the  tailwater.  The 
minimum  water-surface  elevation  at  sta  17+00  was  302.9  and  occurred 
approximately  8  min  after  the  emptying  valve  began  opening.  The  water-surface 
differential  between  these  stations  is  shown  in  Figure  15.  The  maximum  positive 
slope  was  0.56  ft  in  600  ft,  and  the  maximum  negative  slope  was  0.59  ft  in  600  ft 
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A  negative  slope  of  0.59  ft  in  600  ft  converted  to  a  hawser  force  as  previously 
described  gives  an  upstream  longitudinal  hawser  force  of  17.0  tons  for  a  3  x  3 
barge  arrangement  with  the  upstream  end  moored  at  sta  17+00.  This  is  only  a 
slight  reduction  from  the  force  of  21.3  tons  computed  for  the  1.5-min  valve  and 
minimum  excavation  plan. 

Velocity  vectors  and  water-surface  contours  computed  for  the  moderate 
excavation  plan  and  a  1.5-min  valve  speed  at  92, 180,  and  360  sec  after  the  valve 
began  opening  are  shown  in  figure  16.  Velocities  of  up  to  5  ft/sec  were  observed 
in  the  lower  approach  channel,  but  this  magnitude  of  velocity  was  not  as 
widespread  as  with  the  minimum  excavation  plan.  Higher  velocities  (>5  ft/sec) 
occurred  in  the  vicinity  of  shallower  depths,  the  left  half  (looking  downstream)  of 
die  lower  approach  channel  between  sta  22+00  and  26+00.  Here,  the  bed 
elevation  was  the  same  as  with  the  minimum  excavation  plan. 

I. 5-min  valve,  sta  17+00  and  27+00.  Additional  time-histories  were  also 
plotted  with  the  1.5-min  valve  speed  to  evaluate  the  water-surface  slopes  between 
sta  17+00  and  27+00  with  die  moderate  excavation  plan.  The  maximum  water- 
surface  differential  was  found  to  be  -0.99  ft  (figure  15),  which  would  result  in  an 
upstream  longitudinal  hawser  force  of  283  tons  with  the  assumptions  stated 
previously. 

II. 7-min  valve,  sta  17+00  and  23+00.  Time-histories  of  the  water-surface 
elevations  at  sta  17+00  and  23+00  with  the  moderate  excavation  plan  and  the 

1 1.7-min  valve  are  shown  in  figure  14.  The  water  surface  varied  only  slightiy 
from  the  tailwaler  elevation,  and  the  maximum  positive  and  negative  water- 
surface  slopes  shown  in  figure  15  were  0.09  and  0.08  fit,  respectively.  These 
slopes  were  slightiy  less  than  those  with  the  minimum  excavation  plan  and  slow 
valve  and  would  result  in  low  hawser  forces.  A  positive  slope  of  0.09  ft  in  600  ft 
would  result  in  a  downstream  longitudinal  hawser  force  of  2.6  tons. 

ll  .7-min  valve,  sta  17+00  and  27+00.  Time-histories  of  the  water  surface 
were  also  obtained  with  the  1 1.7-min  valve  speed  with  the  moderate  excavation 
plan  In  particular,  die  water-surface  slopes  between  sta  17+00  and  27+00  were 
examined.  The  maximum  water-surface  differential  was  found  to  be  -0. 14  ft 
(figure  15),  which  would  result  in  an  upstream  longitudinal  hawser  force  of 
4.0  tons  with  the  assumptions  stated  previously. 

Comparison  of  results  with  moderate  excavation.  Comparisons  of  the 
results  with  13-  and  11.7-min  valve  times  given  the  moderate  excavation  plan  are 
also  shown  in  figure  15.  The  faster  valve  speed  (13  min)  produced  much  larger 
water-surface  slopes  than  the  slower  (11.7  min)  valve,  and  these  higher  slopes 
would  cause  higher  hawser  forces  for  a  tow  and  barges  moored  in  the  lower  lock 
approach.  The  flow  conditions  in  the  lower  approach  were  favorable  with  die 
1 1.7-min  valve,  but  slower  emptying  times  would  result  from  the  slower  valve. 
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Comparison  of  minimum  and  moderate  excavation  plans 

A  comparison  of  Figures  1 1  and  15  shows  the  difference  between  the  water- 
surface  differentials  for  the  minimum  and  moderate  excavation  plans  for  the 
1.5-min  valve.  The  maximum  positive  and  negative  differentials  are  reduced 
slightly  with  the  moderate  excavation  plan  and  result  from  the  larger  depths  with 
this  plan.  However,  this  reduction  is  not  enough  to  significantly  reduce  the 
hawser  forces  for  a  tow  moored  in  the  lower  approach.  A  very  slight  reduction  in 
the  maximum  and  minimum  water-surface  differentials  was  also  observed  with 
the  1 1 ,7-min  valve  and  moderate  excavation  plan  as  seen  by  comparing 
Figures  1 1  and  15.  The  water-surface  differentials  observed  for  the  minimum  and 
moderate  excavation  plans  with  the  1 1 ,7-min  valve  were  essentially  the  same. 
Results  from  these  simulations  are  provided  in  Table  6. 


Landside  Discharge  Channel  Alternative 

The  Nashville  District  requested  that  WES  evaluate  flow  conditions  in  the 
lower  lock  approach  with  a  landside  channel  that  discharged  the  lock  flow  at  the 
end  of  the  landside  guide  walL  A  plan  view  of  the  landside  discharge  channel  is 
shown  in  Figure  3.  The  following  flow  conditions  were  evaluated  for  this 
discharge  alternative: 
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The  discharge  hydrographs  for  these  emptying  conditions  at  sta  26+20  in  the 
landside  channel  were  furnished  by  the  Nashville  District  and  are  shown  in 
Figures  8  and  9. 

Simulations  were  made  using  3-sec  time-steps.  As  with  the  lateral  discharge 
design,  the  total  simulation  times  were  limited  by  file  total  time  provided  on  the 
inflow  hydrographs.  The  computational  mesh  used  for  all  simulations  with  the 
landside  discharge  alternative  is  shown  in  Figure  6.  The  mesh  represented 
approximately  one  mile  of  topography  downstream  from  file  proposed  lock  and  a 
1,600-ft  width  of  file  river  at  the  discharge  outlet  Downstream  from  the 
discharge  outlet,  the  mesh  represented  approximately  a  1,250-ft  width  of  the 
river. 
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Headwater  el  359,  taifwater  el  302, 
1.5-min  valve,  sta  17+00  and  23+00 


The  first  simulation  with  the  landside  discharge  channel  alternative  was 
performed  with  headwater  el  359,  tailwater  el  302,  and  a  1.5-min  emptying  valve 
operation.  Water-surface  elevations  in  the  lower  lock  approach  were  compared 
to  evaluate  the  water-surface  slopes.  The  locations  chosen  were  similar  to  those 
locations  used  in  the  evaluation  of  the  flow  conditions  with  the  interlaced  lateral 
discharge  system.  The  upstream  location  is  at  sta  17+00,  the  downstream 
location  is  at  sta  23+00,  and  both  are  55  ft  out  from  the  landside  guide  wall.  The 
minimum  excavation  plan  was  used  for  the  topography  in  the  lower  lock 
approach.  Time-histories  of  the  water-surface  elevations  at  sta  17+00  and  23+00 
are  shown  in  Figure  17.  The  water-surface  differential  between  these  stations, 
shown  in  Figure  18,  represents  the  water-surface  slope  between  the  two  locations. 
The  maximum  differential  occurred  160  sec  after  the  hydrograph  began  and  was 
0.69  ft.  A  positive  differential  means  the  upstream  water  surface  is  higher  than 
the  downstream  water  surface,  and  a  negative  differential  indicates  the 
downstream  water  surface  is  higher.  If  one  neglects  the  forces  due  to  drag  and 
assumes  that  the  force  on  the  barges  results  entirely  from  the  water-surface  slope, 
the  longitudinal  hawser  force  computed  for  a  3  x  3  barge  arrangement  with  a  9-ft 
draft  and  a  water-surface  slope  of  0.69  ft  in  600  ft  is  20.1  tons.  The  longitudinal 
hawser  force  computed  for  the  minimum  excavation  plan,  headwater  el  357, 
tailwater  el  304.2,  and  a  1.5-min  valve  operation  with  the  interlaced  lateral 
discharge  alternative  was  21.3  tons. 

Velocity  vectors  and  water-surface  contours  at  90, 180,  and  360  sec  after  the 
valve  began  opening  are  shown  in  Figure  19.  Large  eddies  in  the  river  are  formed 
by  the  high-velocity  flow  discharging  into  the  river  as  shown  in  Figure  19b. 

Large  water-surface  gradients  also  occurring  in  the  vicinity  of  these  eddies  are  not 
a  desirable  hydraulic  feature.  Swirling  high-velocity  flow  in  the  river  could  be  a 
safety  hazard  for  small  craft  as  well  as  cause  navigation  difficulties  for  larger 
tows. 


Headwater  el  359,  tailwater  el  302, 

1.5-min  valve,  sta  15+00  and  27+00 

Flow  conditions  in  the  lower  lock  approach  were  also  evaluated  at  sta  15+00 
and  27+00  to  determine  the  longitudinal  water-surface  slopes  that  a  3  x  6  barge 
arrangement  would  experience  with  the  upstream  end  at  sta  15+00.  Time- 
histories  of  the  water-surface  elevations  at  these  stations  with  headwater  el  359, 
tailwater  el  302,  and  a  15-min  emptying  valve  are  shown  in  Figure  20.  A  time- 
history  of  the  water-surface  differential  between  these  stations  is  shown  in 
Figure  21.  The  maximum  water-surface  differential  was  1.32  ft  and  occurred  2.3 
min  after  the  valve  began  to  open.  For  the  assumptions  already  stated,  this  slope 
results  in  a  downstream  longitudinal  hawser  force  of  37.9  tons. 
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Headwater  el  359,  tailwater  el  314, 1.5-min  valve, 
sta  17+00  and  23+00 

The  next  simulation  was  performed  with  headwater  el  359,  tailwater  el  314, 
and  a  1.5-min  valve  operation-  Time-histories  of  the  water-surface  elevations  at 
the  upstream  and  downstream  stations  with  these  conditions  are  shown  in 
Figure  17.  The  maximum  water-surface  differential  was  0.37  ft  (Figure  18)  and 
occurred  2.6  min  after  the  emptying  valve  began  to  open.  A  water-surface  slope 
of  0.37  ft  per  600  ft  for  a  3  x  3  barge  arrangement  drafted  to  9  ft  with  the 
assumptions  stated  gives  a  hawser  force  of  10.6  tons.  This  is  almost  one-half  the 
force  determined  with  tailwater  el  302.  Greater  depths  associated  with  a  higher 
tailwater  change  the  characteristics  of  the  lower  lock  approach  channel.  These 
increased  depths  result  in  faster  wave  celerity  (gk) 1/2  and  shorter  periods  of  water- 
surface  oscillations  within  the  channel.  Although  the  magnitude  of  the 
oscillations  did  not  differ  significantly  from  those  of  a  1.5-min  valve  time  and 
tailwater  el  302  (Figure  17),  the  differentials  between  the  stations  within  the 
channel  were  considerably  less.  The  water-surface  differential  between  the 
upstream  and  downstream  stations  and  consequently  the  hawser  forces  are  less. 


Headwater  el  359,  tailwater  el  314, 1.5-min  valve, 
sta  15+00  and  27+00 

Time-histories  of  die  water-surface  elevations  at  sta  15+00  and  27+00  with 
headwater  el  359,  tailwater  el  314,  and  a  1.5-min  emptying  valve  are  shown  in 
Figure  20.  A  time-history  of  die  water-surface  differential  between  these  stations 
is  shown  in  Figure  21.  The  maximum  water-surface  differential  occurred  2.6  min 
after  the  valve  began  to  open.  This  differential  of  0.88  ft  results  in  a  downstream 
longitudinal  hawser  force  of 253  tons  using  die  assumptions  mentioned 
previously. 


Headwater  el  359,  tailwater  el  302, 6.0-min  valve, 
sta  17+00  and  23+00 

A  simulation  was  performed  next  with  headwater  el  359,  tailwater  el  302,  and 
a  6-min  emptying  valve  operation.  The  time-histories  of  water-surface  elevation 
at  sta  17+00  and  23+00  are  shown  in  Figure  22.  The  water-surface  differential 
between  the  stations  is  shown  in  Figure  18.  The  maximum  water-surface 
differential  was  -0.1 1  ft  and  occurred  at  1 1.3  min  during  the  emptying  cycle. 

This  differential  gives  a  hawser  force  of  2.6  tons.  The  negative  differential 
indicates  the  downstream  water  surface  was  higher  than  the  upstream  water 
surface. 
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Headwater  el  359,  tailwater  el  302, 6.0-min  valve, 
sta  15+00  and  27+00 

Time-histories  of  the  water-surface  elevations  at  sta  15+00  and  27+00  with 
headwater  el  359,  tailwater  el  302,  and  a  6-min  emptying  valve  are  shown  in 
Figure  23.  The  water-surface  differential  between  these  stations  is  plotted  as  a 
function  of  time  in  Figure  21.  The  maximum  water-surface  differential  was 
0.22  ft  and  occurred  3.6  min  after  the  valve  began  to  open.  For  the  assumptions 
already  stated,  this  slope  results  in  a  downstream  longitudinal  hawser  force  of 
6.3  tons. 


Comparison  of  landside  discharge  channel  alternative  results 

The  results  from  the  three  simulations  are  summarized  in  Table  7.  These 
results  indicate  that  the  slower  valve  speed  has  a  significant  influence  on  the 
water-surface  slopes  in  the  lower  lock  approach  as  one  would  expect  The 
hawser  forces  computed  based  on  longitudinal  water-surface  slope  were  reduced 
by  almost  an  order  of  magnitude  by  slowing  die  valve  opening  time  from  1.5  min 
to  6.0  min.  The  hawser  forces  were  also  reduced  significantly  with  a  12-ft  higher 
tailwater  for  the  1-5-min  emptying  valve  operation.  It  is  interesting  to  note  that 
even  though  die  landside  channel  discharges  downstream  from  a  tow  moored 
along  the  landside  guide  wall,  for  five  of  die  six  simulations,  the  longitudinal 
hawser  forces  are  in  the  downstream  direction.  This  is  a  result  of  the  water- 
surface  oscillations  caused  by  the  long-period  gravity  wave  generated  from  the 
lock  discharging  into  the  river.  Maximum  oscillations  are  present  at  the  vertical 
barrier  of  the  lock  lower  miter  gates. 

Because  die  maximum  drawdown  at  the  lower  miter  gates  occurs  relatively 
early  in  the  emptying  cycle,  there  should  still  be  positive  head  on  the  miter  gates. 
However,  conditions  of  reverse  head  on  the  miter  gates  are  possible  if  the 
oscillations  remain  after  die  lock  chamber  has  emptied.  Successful  design  of  the 
miter  gate  mechanical  components  depends  on  accurate  determinations  of  the 
pressures  on  the  miter  gales. 
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5  Summary  and  Conclusions 


Two-dimensional  (depth-averaged)  simulations  of  the  flow  conditions  in  the 
lower  lock  approach  for  the  proposed  1,200-ft  lock  at  Kentucky  Lock,  Tennessee 
River,  were  performed  to  evaluate  the  interlaced  lateral  and  landside  channel 
discharge  alternatives.  Two  bed  configurations  were  evaluated  with  the 
interlaced  lateral  discharge  alternative  to  determine  the  benefits  of  increased 
excavation  in  die  lower  approach  channel.  The  simulations  with  the  interlaced 
lateral  discharge  system  were  conducted  for  headwater  el  357,  tailwater  el  304.2, 
and  emptying  valve  opening  times  of  13  and  1 1.7  mrn.  This  headwater  and 
tailwater  combination  represents  the  50  percent  duration  conditions,  and  the 
valve  times  represent  fast  and  slow  valve  operations.  The  simulations  with  the 
landside  discharge  channel  were  conducted  with  headwater  el  359,  tailwater 
el  302  and  314,  and  a  13-min  emptying  valve  operation.  An  emptying  valve 
opening  time  of  6  min  was  also  simulated  with  headwater  el  359  and  tailwater 
el  302. 

The  flow  conditions  in  the  lower  approach  were  evaluated  by  comparing 
water-surface  differentials  at  selected  locations.  These  locations  were  in  the 
vicinity  of  the  bow  and  stem  of  a  3  x  3  barge  arrangement  moored  with  the 
upstream  end  at  sta  17+00  and  a  3  x  5  barge  arrangement  with  the  upstream  end 
moored  at  sta  17+00  for  the  interlaced  lateral  discharge  alternative.  With  the 
landside  discharge  alternative,  flow  conditions  were  also  evaluated  at  sta  15+00 
and  27+00,  which  would  be  representative  of  a  3  x  6  barge  arrangement  with  the 
upstream  end  moored  at  sta  15+00.  The  water-surface  slope  between  the  two 
stations  is  considered  to  be  a  major  contributor  to  the  hawser  force  the  barges 
moored  in  this  area  would  experience  during  lock  emptying.  In  this  type  of 
analysis,  the  hawser  forces  the  tow  and  barge  arrangement  will  experience  are 
directly  related  to  the  water-surface  slope  on  which  the  vessel  rests. 

The  simulation  results  with  the  interlaced  lateral  discharge  system  indicate 
that  with  the  minimum  excavation  plan  and  the  fast  (13  min)  valve,  large 
longitudinal  water-surface  slopes  exist  in  the  lower  approach.  A  3  x  3  barge 
arrangement  with  the  upstream  end  moored  at  sta  17+00  would  experience 
hawser  forces  greater  than  20  tons.  A  3  x  5  barge  arrangement  with  the  upstream 
end  moored  at  sta  17+00  would  experience  hawser  forces  greater  than  34  tons. 

The  water-surface  slopes  resulting  from  the  1 1.7-min  valve  and  minimum 
excavation  were  greatly  reduced  from  those  observed  with  the  fast  valve.  A  3  x  3 
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and  a  3  x  5  barge  arrangement  moored  with  the  upstream  end  at  sta  17+00  would 
probably  experience  hawser  forces  less  than  5  tons  with  the  1 1.7-min  emptying 
valve.  Of  course,  slow  valve  operations  result  in  slower  emptying  times,  which 
must  be  considered  as  well  as  the  flow  conditions  inside  the  lock  chamber  and 
culverts  to  determine  the  overall  effects  from  slow  valve  times. 

The  water-surface  slopes  computed  for  the  1.5-  and  1 1.7-min  valve  speeds 
with  the  interlaced  lateral  discharge  alternative  and  the  moderate  excavation  plan 
were  slightly  less  than  those  computed  with  the  minimum  excavation  plan,  as  one 
would  expect  given  the  same  boundary  conditions.  The  surge  produced  by  the 
lock  discharge  is  essentially  independent  of  the  bed  elevation,  although  the 
greater  depths  resulting  from  the  moderate  excavation  would  coincide  with 
slower  velocities.  The  seiching  frequency  within  the  approach  channel  is  only 
slightly  higher  with  the  5-ft  increase  in  depth.  The  gravity  wave  celerity  is 
increased  by  a  factor  of  the  square  root  of  the  moderate  excavation  depth  (20.2  ft) 
to  minimum  excavation  depth  (15.2  ft)  ratio.  That  is,  the  moderate  excavation 
plan  increases  the  celerity  by  a  factor  of  (20-2/15.2)  m  or  1.15.  This  increase  in 
celerity  results  in  a  reduction  of  the  oscillation  period  (increased  frequency)  of  the 
lower  approach. 

Relatively  large  slopes  still  occurred  with  the  moderate  excavation  plan  with 
the  fast  valve,  and  much  smaller  slopes  occurred  with  the  slow  valve.  The 
hawser  forces  computed  for  die  moderate  excavation  plan  were  less  than  those 
computed  with  die  minimum  excavation  plan,  but  not  significantly  less.  The 
additional  5  ft  of  excavation  in  die  moderate  excavation  plan  does  not  provide 
enough  reduction  in  hawser  force  to  consider  this  plan.  The  differences  observed 
between  the  two  excavation  plans  indicate  that  a  certain  amount  of  excavation  is 
probably  beneficial.  This  amount  should  probably  be  based  on  the  amount 
needed  for  tow  and  barge  maneuverability  in  the  lower  approach  including  vessel 
squat,  rather  than  on  an  amount  needed  to  reduce  hawser  forces.  Franco  (1976) 
points  out  that  if  the  channel  is  not  deep  enough,  barges  can  hit  the  bottom  of  the 
canal  during  the  wave  troughs. 

The  results  from  the  simulations  with  die  landside  discharge  channel 
alternative  support  the  results  obtained  with  the  interlaced  lateral  discharge 
alternative.  A  slower  emptying  valve  reduces  the  water-surface  slopes  in  the 
lower  lock  approach  for  die  same  headwater  and  tailwater  combination.  A  higher 
tailwater  elevation  helps  reduce  the  hawser  forces,  but  die  slower  valve  is  where 
die  most  reduction  in  hawser  force  can  be  achieved.  The  results  computed  from 
the  interlaced  lateral  discharge  alternative  and  the  landside  discharge  channel 
alternative  are  not  directly  comparable  due  to  the  different  headwater  and 
tailwater  combinations  requested  by  the  Nashville  District  However,  for  the 
same  headwater  and  tailwater  conditions,  the  interlaced  lateral  system  would 
probably  result  in  slightly  higher  hawser  forces  for  tow  and  barge  arrangements 
moored  along  the  landside  guide  waLL  A  3  x  6  barge  arrangement  with  the 
upstream  end  moored  at  sta  17+00  would  probably  experience  hawser  forces  over 
40  tons  with  the  interlaced  laterals  compared  with  forces  greater  than  37  tons  for 
a  3  x  6  with  the  upstream  end  moored  at  sta  15+00  with  the  landside  discharge 
channel  alternative.  This  difference  is  relatively  small,  and  does  not  justify 
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construction  of  the  higher  cost  landside  channel.  The  most  effective  lock  outlet 
from  a  hydraulic  engineer’s  viewpoint  is  one  that  is  hydraulically  efficient  yet 
provides  good  energy  dissipation  and  has  minimal  impact  on  navigation. 
Unfortunately,  the  economics  involved  for  this  type  system  usually  eliminates  it 
as  an  alternative.  The  performance  of  the  outlet  must  then  be  adjusted  to  satisfy 
the  economics.  Discharging  a  lock  immediately  downstream  from  the  lower 
miter  gates  with  an  interlaced  lateral  is  generally  less  expensive,  but  requires 
extreme  caution  when  operating  the  emptying  valves.  Slower  valve  operations 
and  lock  emptying  times  will  be  required  for  this  type  discharge  system  versus 
one  that  discharges  away  from  the  lower  approach. 
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Figure  1 .  Vicinity  map 


Figure  2.  Interlaced  lateral  discharge  system 
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Figure  3.  Landside  discharge  channel 
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Figure  4.  Plan  view  of  lower  lock  approach,  interlaced  lateral  discharge  alternative  and  mini 


a.  Overall  view  of  mesh 

Figure  5.  Mesh  used  for  hydrodynamic  computations  with  interlaced  lateral  discharge  alternative  (continued) 


b.  Close-up  of  lower  approach 
Figure  5.  (Concluded) 
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Figure  6.  Computational  mesh  used  for  simulations  with  the  landside  discharge  channel 


Figure  7.  Emptying  valve  hydrographs,  headwater  el  357,  tailwater  el  304.2,  interlaced  lateral  discharge 
alternative  boundary  conditions 


Figure  8.  Emptying  valve  hydrograph,  sta  26+20,  1 .5-min  valve  time,  headwater  el  359,  landside 
discharge  channel  alternative 


Figure  9.  Emptying  valve  hydrograph,  sta  26+20,  6-min  valve,  headwater  el  359,  landside  discharge 
channel  alternative 
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Figure  10.  Water-surface  elevations,  tailwater  el  304.2,  minimum  excavation  plan,  interlaced  lateral  discharge  alternative 


Figure  1 1 .  Water-surface  differentials,  lower  pool  el  304.2,  minimum  excavation  plan,  interlaced  lateral  discharge  alternative 


a.  Conditions  at  92  sec 


Figure  12.  Velocity  vectors  and  water-surface  contours  during  lock  emptying,  lower  pool  el  304.2, 
1 ,5-min  valve,  minimum  excavation  plan,  interlaced  lateral  discharge  alternative 
(Sheet  1  of  3) 


b.  Conditions  at  1 80  sec 


Figure  1 2.  (Sheet  2  of  3) 


c.  Conditions  at  360  sec 


Figure  1 2.  (Sheet  3  of  3) 
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13.  Moderate  excavation  plan,  interlaced  lateral  discharge  alternative 
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Figure  14.  Water-surface  elevations,  1 ,5-min  valve,  lower  pool  el  304.2,  moderate  excavation  plan,  interlaced  lateral  discharge  alternative 


a.  Conditions  at  92  sec 


Figure  16.  Velocity  vectors  and  water-surface  contours  during  lock  emptying,  lower  pool  el  304.2, 
1 ,5-min  valve,  moderate  excavation  plan,  interlaced  lateral  discharge  alternative 
(Sheet  1  of  3) 


b.  Conditions  at  1 80  sec 


Figure  1 6.  (Sheet  2  of  3) 


c.  Conditions  at  360  sec 


Figure  1 6.  (Sheet  3  of  3) 
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Figure  17.  Water-surface  elevations,  sta  17+00  and  23+00,  headwater  el  359, 1 .5-min  valve,  landside  discharge  channel  alternative 


STA  17+00 


Figure  18.  Water-surface  differentials  between  sta  17+00  and  23+00,  landside  discharge  channel  alternative 


a.  Conditions  at  90  sec 

Figure  19.  Velocity  vectors  and  water-surface  contours  during  lock  emptying,  lower  pool  el  302, 
1 ,5-min  valve,  landside  discharge  channel  alternative  (Sheet  1  of  3) 


b.  Conditions  at  180  sec 


Figure  1 9.  (Sheet  2  of  3) 


c.  Conditions  at  360  sec 


Figure  1 9.  (Sheet  3  of  3) 


Figure  20.  Water-surface  elevations,  headwater  el  359,  sta  1 5+00  and  27+00, 1 ,5-min  valve,  landside  discharge  channel  alternative 


e  22.  Water-surface  elevations,  sta  17+00  and  23+00,  headwater  el  359,  tailwater  el  302,  6-min  valve,  landside  discharge  channel 


e  23.  Water-surface  elevations,  sta  15+00  and  27+00,  headwater  el  359,  tailwater  el  302,  6-min  valve,  landside  discharge  channel 
alternative 
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Table  1  (Continued! 


Time,  sec 

444 

446 

448 

450 

452 

454 

456 

I458 

460 

462 

1  464 

466 

|  468 

l470 

I472 

474 

476 

478 

I  480 

482 

484 

486 

488 

490 

492 

494 

496 

498 

500 

502 

. 

b — 

Total  Pischai 


9300.7 


9266.8 


9232.9 


9198.9 


9164.9 


9131 


9097.1 


9063.1 


9029.2 


8995.4 


8961.5 


8927.7 


8893.9 


8860.1 


8826.3 


8792.5 

8758.7 

8724.9 


8691 


8657.1 

8623.1 


8589.1 


8555.1 


8521 


8487 


8452.9 


8418.8 


8384.7 


8350.6 


8316.5 


8282.5 


506 


8248.4 
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Table  1  (Continued 


Time,  sec 

572 

574 

576 

578 

580 

582 

584 

586 

5 88 

590 

592 

594 

596 

598 

600 

602 

tt  604 

|  606 

608 

610 

612 

614 

616 

618 

620 

622 

624 

626 

628 

630 

I632 

B 

Total  Dischat 

7128.7 

7094.7 

7060.6 

7026.6 

6992.5 

6958.4 

6924.4 

6890.3 

6856.3 

6822.3 

6788.3 

6754.3 

6720.3 

6686.3 

6652.3 

6618.3 

6584.4 

6550.4 

6516.4 

6482.4 

6448.5 

6414.5 

6380.6 

6346.6 

6312.7 

6278.8 

6244.8 

6210.9 

6177.1 

6143.2 

6109.3 


634 


6075.4 


Table  1  (Continued 


Time,  sec 

Total  Discharge,  cfs 

636 

6041.6 

638 

6007.7 

640 

5973.8 

642 

5939.9 

644 

5906 

646 

5872.1 

648 

5838.2 

650 

5804.2 

652 

5770.3 

654 

5736.3 

656 

5702.3 

658 

5668.3 

660 

5634.2 

662 

5600.2 

664 

5566.1 

666 

5532.1 

668 

5498 

670 

5464 

672 

5429.9 

674 

5395.9 

676 

5361.8 

678 

5327.8 

680 

5293.8 

682 

5259.7 

684 

5225.7 

686 

5191.7 

690 

5157.7 

692 

5123.7 

694 

5089.6 

696 

5021.6 

698 

4987.6 

700 

4953.6 

i 
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Table  2  (Continued] 


ji 


Time,  sec 

Total  Discharge,  cfs 

252 

4442.6 

254 

4507.5 

256 

4575.4 

258 

4646.7 

260 

4704.3 

262 

4747.6 

264 

4792 

266 

4840.2 

268 

4890.7 

270 

4942.2 

272 

4993.8 

274 

5045.4 

_ — - 

276 

5097.1 

278 

5149.7 

280 

5203.7 

282 

5259.7 

284 

5318 

286 

5378.9 

28 8 

5442.5 

290 

5508.7 

292 

5574.6 

294 

5627.5 

296 

5671.4 

298 

5717.7 

300 

5767.2 

302 

5818.5 

304 

5870.6 

|  306 

5922.7 

308 

5974.9 

310 

6027.3 

312 

6080.7 

314 

6135.3 

[Sheet  5  of  20] 


Table  2  (Continued' 

Time,  sec 

508 

510 

512 

514 

516 

518 

520 

522 

524 

526 

528 

530 

532 

534 

536 

538 

540 

542 

544 

546 

548 

550 

552 

554 

556 

558 

560 

562 

564 

566 

568 

10747 

10771.8 

10796.5 

10820.4 

10843.3 

10865.6 

10887.8 

10910.3 

10933.3 
10952 

10964.1 

10973.1 
10982 

10992.3 

11004.2 

11016.7 

11028.7 

11039.3 

11048.4 

11056.4 
11064 

11071.8 

11079.8 

11085.7 

11088.7 

11090.2 

11091.6 

11093.4 

11095.7 

11098.1 

11100.2 


570 


11101.6 


Table  2  (Continued 


Time,  sec 

Total  Discharge,  cfs 

572 

11102.3 

574 

11102.4 

576 

11102.4 

578 

11102.4 

580 

11102 

582 

11099.3 

584 

11094.1 

586 

11087.8 

588 

11081.4 

590 

11075.6 

592 

1 1070.2 

594 

11064.9 

596 

11059.1 

598 

11052.5 

600 

11044.9 

602 

11036.5 

604 

11027.7 

606 

11018.9 

608 

11009.6 

610 

10997.9 

612 

10983.5 

614 

10967.9 

616 

10952.1 

618 

10936.9 

620 

10922.1 

622 

10907.4 

624 

10892.1 

626 

10876 

628 

10859 

630 

10841.2 

632 

10823 

634 

10804.8 
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698 


10069 


Table  2  (Continued] 

Time,  sec 

764 

766 

768 

770 

772 

774 

776 

778 

780 

782 

784 

786 

788 

790 

792 

794 

796 

798 

800 

802 

804 

806 

808 

810 

812 

814 

816 

818 

820 

822 

824 

1 


Total  Dischai 

8966.8 

8932.9 

8899 

8865 

8831 

8797 

8763 

8729 

8695 

8661 

8627 

8593.1 

8559.1 

8525.1 

8491.2 

8457.3 

8423.4 

8389.5 

8355.6 

8321.7 

8287.8 

8253.8 

8219.9 

8185.9 

8152 

8118 

8084 

8050 

8016 

7982 

7984 

826 


7914.1 


828 

7880.1 

830 

7846.2 

832 

7812.3 

834 

7778.3 

836 

7744.4 

83 8 

7710.5 

840 

7676.6 

842 

7642.6 

844 

7608.7 

846 

7574.7 

848 

7540.8 

850 

7506.8 

852 

7472.8 

854  7438.8 

856  7404.8 

858 
860 
862 
864 
866 
868 
870 
872 
874 
876 
878 
880 
882 
884 
686 
8 88 
890 


1 


7370.8 

7336.8 

7302.8 
7268 
7234 

7200.9 
7167 
7133 

7099.1 

7065.1 

7031.2 

6997.3 

6963.4 

6929.4 

6895.5 

6861.6 
6827.6 
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Table  3 

Hydrograph  and  Stages  at  Sta  26+20, 1.5-min  Emptying  Valve,  Headwater  EL  359, 
Tailwater  El  302,  Landside  Discharge  Channel  Alternative 


Time,  sec 


Total  Discharge,  cfs 


Table  4 

Hydrograph  and  Stages  at  Sta  26+20, 1.5-min  Emptying  Valve,  Headwater 
El  359,  Tailwater  EL  314,landside  Discharge  Channel  Alternative _ 


Time,  sec 


Total  Dischar 


Stage,  feet 


Table  4  (Continued 


Time,  sec 


Total  Discharge,  cfs 


Table  4  (Continued 


Time,  sec 

Total  Discharge,  cfs 

Stage,  feet 

184 

17587.51 

312.7 

186 

17306.71 

312.7 

188 

17026.22 

312.8 

190 

16747.62 

312.8 

192 

16471.21 

312.8 

194 

16196.01 

312.9 

196 

15919.79 

312.9 

198 

15639.69 

313.0 

200 

15353.54 

313.0 

202 

15060.67 

313.1 

I204 

14761.4 

313.1 

206 

14455.62 

313.1 

208 

14142.28 

313.2 

210 

13820.17 

313.2 

212 

13488.79 

313.3 

214 

13148.76 

313.3 

216 

12802.1 

313.3 

218 

12452.58 

313.4 

220 

12105.58 

313.4 

222 

11767.31 

313.5 

224 

11443.6 

313.5 

226 

11138.6 

313.5 

228 

10853.67 

313.5 

230 

10587.02 

313.6 

232 

10334.44 

313.6 

234 

10090.84 

313.6 

236 

9852.238 

313.6 

238 

9616.955 

313.6 

240 

9385.33 

313.7 

242 

9158.299 

313.7 

244  1 

8936.149 

313.7 

Table  4  (Continued 


Time,  sec 


Total  Discharge,  cfs 


Stage,  feet 


5851.655 


Table  4  (Concluded 


Time,  sec 


Total  Discharge,  cfs 


Stage,  feet 


1542.158 


Table  5 

Hydrograph  and 
Tailwater  El  302, 


Stages  at  Sta  26+20,  6-min  Emptying  Valve,  Headwater  El  359, 
Lanside  Discharge  Channel  Alternative 


Time,  sec 


Total  Dischar 


Table  5  (Continued 


Time,  sec 


Total  Discharge,  cfs 


3071.019 


Time,  sec _  Total  Discharge,  cfs  Stage,  feet 

122  3869.528  301.8 


124 

3995.885 

301.8 

126 

4138.977 

301.7 

128 

4301.425 

301.7 

130 

4474.812 

301.7 

132 

4650.235 

301.7 

134 

4824.212 

301.7 

136 

4993.595 

301.6 

138 

5149.691 

301 .6 

140 

5282.091 

301.6 

142 

5388.838 

301.6 

144 

5481.111 

301.5 

146 

5576.836 

301.5 

148 

5688.853 

301.5 

150 

5817.425 

301.5 

I  152 

5952.904 

301.5 

|  154 

6085.781 

301.4 

1  156 

6213.664 

301.4 

1 158 

6338.299 

301.4 

1 160 

6459.741 

301.4 

162 

6578.189 

301.3 

164 

6699.549 

301.3 

166 

6833.009 

301.3 

168 

6981.971 

301.2 

170 

7140.614 

301.2 

172 

7299.106 

301.2 

174 

7448.829 

301.1  1 

176 

7583.698 

301.1 

178 

7702.162 

301.0 

180 

7810.455 

301.0 

182 


7920.761 


301.0 


Table  5  (Continued 


Time,  sec 


Total  Discharge,  cfs 


Stage,  feet 
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496 


12616.35 


298.3 


498 

12569.89 

298.4 

500 

12523.78 

298.4 

502 

12477.93 

298.5 

504 

12432.43 

298.5 

506 

12387.24 

298.6 

508 

12342.36 

298.6 

510 

12297.71 

298.7 

512 

12253.32 

298.7 

514 

12209.18 

298.8 

516 

12165.3 

298.8 

518 

12121.66 

298.8 

520 

12078.25 

298.9 

522 

12034.92 

298.9 

524 

11991.5 

299.0 

526 

11947.81 

299.0 

528 

11903.82 

299.0 

530 

11859.57 

299.1 

532 

11815 

299.1 

534 

11770.58 

299.1 

536 

11725.81 

299.2 

538 

11680.76 

299.2 

540 

11635.4 

299.2 

542 

11589.76 

299.3 

544 

11543.84 

299.3 

546 

11497.67 

299.3 

548 

11451.29 

299.4 

550 

11404.69 

299.4 

552 

11357.83 

299.4 

554 

11310.7 

299.5 

Table  5  (Continued 


Time,  sec 


Total  Discharge,  cfs 


9876.649 


Table  5  (Continued 


Table  5  (Continued 


Time,  sec 


Total  Discha 


7195.253 


7153.982 


7114.795 


6988.235 


6963.317 


6939.618 


6632.007 


6553.534 


6297.275 
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Time,  sec 

Total  Disharge,  cfs 

Stage,  feet 

804 

6117.593 

301.4 

806 

6073.208 

301.4 

808 

6029.223 

301.4 

810 

5985.633 

301.5 

812 

5942.399 

301.5 

814 

5899.497 

301.5 

816 

5856.945 

301.5 

818 

5814.802 

301.5 

820 

5773.129 

301.5 

822 

5731.948 

301.5 

I824 

5691.233 

301.5 

826 

5650.921 

301.5 

828 

5610.942 

301.5 

830 

5571.244 

301.5 

832 

5531.801 

301.5 

834 

5492.605 

301.5 

836 

5453.64 

301.6 

838 

5414.84 8 

301.6 

840 

5376.108 

301.6 

842 

5337.231 

301.6 

844 

5297.965 

301.6 

846 

5258.023 

301.6 

84 8 

5217.116 

301.6 

850 

5174.98 

301.6 

852 

5131.416 

301.6 

854 

5086.313 

301.6 

856 

.  -  ...  - . - .. 

5039.68 

301.6 

858 

4991.668 

301.6 

860 

4942.571 

301.6 

862 

4892.813 

301.6 

864 

4842.904 

301.7 
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Table  6 

Computed  Hawser  Forces,  Headwater  El  357,  Tailwater  El  304.2,  Interlaced  Lateral 
Discharge  Alternative 


Excavation  Plan 

Valve  Opening 
Time,  min 

Water-Surface  Differential,  ft 

Computed  Longitudinal 

Hawser  Force,  tons 

Sta  17+00  and 
23+00 

Sta  17+00  and 
27+00 

Sta  17+00  and 
23+00 

Sta  17+00  and 
27+00 

Minimum 

1.5 

-0.74 

-1.19 

21.3 

34.2 

11.7 

0.10 

-0.16 

2.9 

4.6 

Moderate 

1.5 

-0.59 

-0.99 

17.0 

28.5 

11.7 

0.09 

-0.14 

2.6 

4.0 

1  Table  7  1 

1  Computed  Hawser  Forces,  Landside  Discharge  Channel  Alternative  I 

Headwater  El 

Tailwater  El 

Valve  Opening 
Time,  min 

Water-Surface  Differential,  ft 

Computed  Longitudinal  I 

Hawser  Force,  tons  1 

Sta  17+00  and 
23+00 

Sta  15+00  and 
27+00 

Sta  17+00  and 
23+00 

Sta  15+00  and  I 
27+00  fl 

|  359 

302  j 

1.5 

0.70 

1.32 

20.1 

37.9  | 

|  359 

314 

1.5 

0.37 

0.88 

10.6 

25.3  H 

||  359 

302 

6.0 

-0.09 

0.22 

2.6 

6.3  i 
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arrangement  would  experience  hawser  forces  greater  than  20  tons,  and  a  3  x  5  barge  arrangement  would  experience 
hawser  forces  grater  than  34  tons.  The  water-surface  slopes  resulting  from  the  1 1.7-min  valve  and  minimum 
excavation  were  greatly  reduced  from  those  observed  with  the  fast  valve.  A  3  x  3  and  a  3  x  5  barge  arrangement 
would  probably  experience  hawser  forces  less  than  5  tons  with  the  1 1.7-min  emptying  valve. 

Evaluation  of  the  interlaced  lateral  design  with  the  moderate  excavation  plan  showed  that  the  water-surface 
slopes  computed  for  the  1.5-  and  1 1.7-min  valve  were  slightly  less  than  those  computed  with  the  minimum 
excavation  plan  as  one  would  expect  given  the  same  boundary  conditions.  The  surge  produced  by  the  lock 
discharge  is  essentially  independent  of  the  bed  elevation,  although  the  greater  depths  resulting  form  the  moderate 
excavation  would  coincide  with  slower  velocities. 

The  results  from  the  simulations  with  the  landside  discharge  channel  alternative  support  the  results  obtained  with 
the  interlaced  lateral  discharge  alternative.  A  slower  emptying  valve  reduces  die  water-surface  slopes  in  the  lower 
lock  approach  for  the  same  headwater  and  tailwater  combination.  A  higher  tailwater  elevation  helps  reduce  the 
hawser  forces,  but  the  slower  valve  is  where  the  most  reduction  in  hawser  force  can  be  achieved.  Discharging  a  lock 
immediately  downstream  from  the  lower  miter  gates  with  an  interlaced  lateral  is  generally  less  expensive,  but 
requires  extreme  caution  when  operating  the  emptying  valves.  Slower  valve  operations  and  lock  emptying  times 
will  be  required  for  this  type  discharge  system  versus  one  that  discharges  away  from  the  lower  approach. 


